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Abstract 

In this paper a four-dimensional hyperchaotic system with only one equilibrium is 
considered and its double Hopf bifurcations are investigated. The general post-bifurcation 
and stability analysis are carried out using the normal form of the system obtained via 
the method of multiple scales. The dynamics of the orbits predicted through the normal 
form comprises possible regimes of periodic solutions, two-period tori, and three-period 
tori in parameter space. 

Moreover, we show how the hyperchaotic synchronization of this system can be re- 
alized via an adaptive control scheme. Numerical simulations are included to show the 
effectiveness of the designed control. 

1 Introduction 

In the last decades study of hyperchaos has received great attention due to its important role 
in nonlinear science \19\ [131 122] • Several hyperchaotic systems with more than one positive 
Lyapunov exponent have been constructed via control techniques. 

In this paper we consider a four dimensional dynamical system firstly introduced in [9j. This 
system was constructed by adding a linear feedback controller to the second equation of the 
chaotic Chen system [5] and an additional new state equation. The rich dynamics of this 
modified driven Chen system has been observed in [10] via computer simulations showing 
both chaotic and hyperchaotic attractors and period-doubling bifurcations. 

The aim of this paper is, from one hand, to understand better the onset of the hyperchaos 
by the investigation of double Hopf bifurcation of the system and the analysis of the post- 
bifurcation dynamics via the normal form theory [T3], [15]. From the other hand we want 
to realize both the chaotic and the hyperchaotic synchronization of the system designing a 
global adaptive controller. 

We will use in Section 3 a perturbation technique, essentially based on multiple scales 
method, to compute the post-double Hopf normal form [231 1211 EH]. This method does not 
need the application of the Center Manifold Theory and can be employed with the aid of 
symbolic computation [261 12Z] • 
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Numerical simulations are performed to corroborate the predictions from the normal form, 
revealing the existence of stable periodic and toroidal attractors in the post-supercritical-Hopf 
cases. 

Moreover, since the interest into synchronization of hyperchaotic systems has been always 
increasing, in particular due to its applications in secure communications (in fact the presence 
of more than one Lyapunov exponent generates more complex dynamics and improves the 
security), we propose a scheme to synchronize the modified hyperchaotic Chen system via 
adaptive control. Recently, many methods and techniques have been developed to realize 
synchronization and control of chaotic systems, see references [T71 [UJ 121 El 03 [29] and [3] 
for a review. In the same spirit of [16], in Section 4 we design an adaptive control law and 
an update rule for uncertain parameters based on Lyapunov stability theory. 

Numerical simulations are presented to verify the effectiveness of the proposed synchroniza- 
tion method both in chaotic and hyperchaotic regime. 



2 Linear analysis and double Hopf bifurcations 

Let us consider the following modified Chen system, firstly obtained in [9]: 

x'i = a(x 2 - xi) 

X2 = —dx\ — X1X3 + CX2 — X4 
X 3 = X\X2 - 6X3 

x\ = xi + k 

where a, b, c, d and k are real constants. 
Once introduced the following coordinates transformations: 



(2.1) 



xi + k, y = yi + k, 



X3 



b ; 



w 



X4 — k [ d — c + 



which translate to the origin the only equilibrium E: 



E = 



k 2 
b 



k 2 



(2.2) 



the system (|2,ip becomes: 



y -- 

z - 
w 



a(y - x) 

—dx — xz + cy — w 
xy — bz — k(x + y) 



k 2 , 
— x + kz 
b 



(2.3) 



By choosing the values of parameters a = 36, b = 3, c = 28 and d = —1Q the origin is unstable 
and W = —a — b + c < 0. Therefore the system is dissipative and the trajectories converge 
to an attractor which is hyperchaotic, as shown in Fig. [TJ 

Several bifurcation routes to hyperchaos from periodic, quasi-periodic and chaotic orbits are 
observed in [TU] via computer simulations. Since there exists only one equilibrium point, the 
possibility of the fixed point undergoing transcritical, pitchfork or saddle-node bifurcations 
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Figure 1: The parameters are chosen as a = 36, b = 3, c = 28, d = — 16 and k = 1. (a) The 
hyperchaotic attractor. (b) The dynamics of the Lyapunov exponents: the system admits 
two positive Lyapunov exponents. 



(all of which involve creation, annihilation or exchange of stability of a pair of equilibria) is 
precluded. Thus the origin can be only destabilized via a Hopf or double Hopf bifurcation. 

Applying the Routh-Hurwitz criterion one gets the following conditions for the origin being 
a stable equilibrium: 

03 > 0, ao > 0, 0302 — a\ > 0, 03(0201 — 0300) — a\ > 0, (2.4) 

where: 



03 = a + b — c 
ak 2 

d2 = ab H — ac + ad — bc+ k 2 

b 

ai = a(-bc + bd + 3k 2 + 1) 
ao = ab 

are the coefficients of the characteristic polynomial p(X) = A 4 + a3A 3 + a2A 2 + aiA + ao- 
Since the post-bifurcation dynamics following regular Hopf bifurcation in this modified Chen 
system has been recently studied in [1], we are interested into double Hopf bifurcations. At 
the double Hopf bifurcation the four eigenvalues of the linearized system must be two pairs 
of purely imaginary complex conjugates Ai,2 = ±^171^3,4 = ±zo"2 and the characteristic 
equation must have the form (A 2 + <r 2 )(A 2 + a^) = 0, which leads to the following conditions 
for a double Hopf bifurcation to occur: 



n c = c = a + b, 

1 



03 

a 2 = a'i + <r% > 0, d = d* = a + b- ^(3k 2 + 1), 

ai = ' 2 , 2 o + 6 3 + 26Vo6 ^ 

ao = aM>0, k> _ , 



«2"4a >0. aua2£ 
both for a, b > and b > 2a or a, b < and b < 2a. 
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3 Normal form and general post-bifurcation dynamics 

To compute the normal form of the modified Chen system (|2.3|) we will use a perturbation 
technique based on the multiple scales method [151 E21 EZ] . 
Near the double Hopf bifurcation different time scales Tj = eH, j = 0,1,2,... can be 
distinguished and therefore the time derivative decouples as follows: 

d d d 2 ® 3 & 

Jt = dr +£ dT 1 +£ W 2 +£ W z '~ { ' 



Let us write the solutions of the original system (|2.3p and the bifurcation parameters as 
nonlinear expansions in e: 



x = exi + e 2 x 2 + e 3 x 3 + 0(e 4 ), y = eyt + e 2 y 2 + e 3 y 3 + 0(e A ), 
z = ezi + e 2 z 2 + e 3 z 3 + 0(e 4 ), w = ew\ + e 2 w 2 + e 3 w 3 + 0(e 4 ), (3.2) 
c = c* + eci + e 2 c 2 + e 3 c 3 + 0(e 4 ), d = d* + edi + e 2 d 2 + e 3 ^3 + 0(e A ), 

where c* , d* are the bifurcation values at the double Hopf singularity as calculated in (|2.5p . 
All the expansion coefficients Xj, y$, z% and i = 1, 2, 3, . . . depend on the time scales Tj, j = 
0,1,2,.... 

Substituting all the above expansions ()3. ![) - ()3.2[) into the system (j2.3[) and collecting the 
terms at each order in e, one gets a sequence of differential systems for the Xi,yi,Zi and 
Wi, i = 1,2,3, ... : 



0(e) 
0(e 2 ) 

0(e 3 ) 

where £ is the following linear operator: 



Cw\ = 0, 
Cw 2 = r 2 , 
Cw 3 = r 3 , 



(3.3) 
(3.4) 
(3.5) 



C 



1 /l 5 4 / 6-c* 
+ 11 + 



A; \ a 5T 4 



+ 



+ (l-6c* + &cf +3^)— + 6 

(7iQ 



+ ~r + b ~ c +d -^Fi 

a b J dTfi 



(3.6) 



and the source terms T 2 and T 3 contains the nonlinear terms. The equations for Xi,yi and z% 
can be easily obtained only in terms of Wi, but here we skip all the calculation details. 
The solution of the linear homogeneous problem (|3.3p is given by: 



u/i = £ K(T fc )e^ To + a j {T h )e-*i T ») , 



(3.7) 



where the fields ctj, depending on the time scales T^, k = 1, 2, . . . , lie on the center manifolds 
and otj represent their complex conjugate fields. The real numbers aj are the imaginary 
parts of the pairs of purely imaginary complex eigenvalues calculated at the double Hopf 
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singularity. Suppressing the secular terms appearing in r 2 only by imposing T\ = and 
c x , d\ =0 one gets the solutions w 2 ■ 

The solvability condition at 0(e 3 ) leads to the following two coupled equations for the fields 
oi\ and a 2 : 

= L i - j) a j + M (j V„ + N^a x a 2 a h (3.8) 

where the coefficients are linearly dependent on the second order deviation c 2 and d 2 
from the bifurcation values. 
Using the polar coordinates aj = Pje l ^j , a>j = pje~ l ^j , with pj and dj depending on T2 and 
separating the real and the imaginary parts, one obtains the following normal form: 

d P± _ „. ( T (J)^ _i_ t _l /i/rO') „2 _l i\r(i) ^2 

r 2 



= Pi ^c 2 + L^d 2 + M^ + iV 1 ^ + 3)J ) (3-9) 



(i) r(i) fl/f (i) „„j Ar(i) 



— - L%c 2 + L%d 2 + M^p] + N? V(- i+3) - (3-10) 



The explicit expression of the coefficients L^ 1 , L y 2 \' , ! and A^ u; , Z = 1,2 in terms of the 
parameters of the original system (|2.3p can be found in [6]. 
The obtained normal form will be an analytical approximation of the periodic orbits at the 
double Hopf singularity. 
The analysis of this normal form closely parallels earlier work by Yu and co-workers, in 
particular [251 E] • 

The stationary solutions of the equations (|3.9p - (|3.10p are the following: 



O : Pl =p 2 = 0, (3.11) 

Hi : p\ = X —r (l$c 2 + L$d 2 ), p 2 = 0, (3.12) 

tfi = I$C2 + L^d 2 + M^pl, 

H 2 : pi = 0, p\ = l -- (l^c 2 + L ( $d 2 ) , (3.13) 



L { £ > c 2 + L^d 2 + M( 2) pl 



2 -Mf > (zffg + 4^) + iVf (iff eg + 4^) 

Pi = _ -I'll - J-X~\ .JV\ 

S (3.14) 



Pi = M^M^-ivWATf) (3 ' 15) 

1?1 = L^l C 2 + 4^2 + M 2 (1) p? + iV^^l, 

tf 2 = 4? c 2 + 4?^ + M ( 2 2) pl + N^pl 

Performing a standard linear analysis of the system (|3.9|) - (|3.1U|) . it follows that the initial 
equilibrium state corresponding to O, is stable if: 
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4? ca + L ( (ld 2 < 0, for j = 1, 2. (3.16) 

The conditions for the existence and stability of the stationary point H\ (and consequently 
of a periodic solution for the system (|2.3p ) are: 

M x (1) < 0, (3.17) 
L { ^c 2 + L^d 2 > 0, (3.18) 

ifl C 2 + 1$ d 2 ~ ~~JTT (4l C 2 + L S <fe) < 0. (3. 19) 

Therefore along the critical line l\\ 

h : 4l c 2 + L^da = 0, 4? c 2 + 4.2 <fc < (3.20) 

the initial equilibrium bifurcates into a family of limit cycle (approximated by the periodic 
solution Hi) as follows from (|3,16p and (|3.18|) . Analogously, H 2 exists stable when: 

M[ 2) < 0, (3.21) 
r(i) 



L { $c 2 + L^d 2 > 0, (3.22) 



4i c 2 + 1$ d 2 - ^- (i® c 2 + 4? d 2 )<0, (3.23) 
M\ ' v ' 

and along the critical line l 2 : 

h ■ 4i c 2 + L^cfc = 0, 4V c 2 + L^cfc < (3.24) 

the initial equilibrium bifurcates into an other family of limit cycles given by the periodic 
solution H 2 . 

Finally, the conditions for the existence and stability of the point T are the following: 

(l^c 2 + L { $d 2 ) - (L<8c2 + L^d 2 ) > (3.25) 



N[ 2) (lQcz + L^dz) - M : (1) (4i Pa + 42*) > ( 3 - 26 ) 

Af^M^ - N^NP > (3.27) 
(L$c 2 + L$d 2 ) (iVf)-M« 

+Af« (LS?c 2 + Lgcfa) (iV« - M?>) < 0. (3.28) 
From the conditions (|3.17| - (|3.19p and (|3,26p it follows that along the line £3: 

h ■■ - ^5y4i J c 2 + ^4? " ~U) L ^J d2 = °' and ~ (M) hold (3.29) 
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the periodic solution Hi bifurcates with frequency $2 into a quasi-periodic solution (with a 
secondary Hopf bifurcation) and a 2-D torus arises. Analogously, taking into account ()3.21|) - 
(|3.23p and (|3.25p . one individuates the line I4: 



k ■ 4i W4? c 2 + iff - -^4? )d2 = 0, and ([321]) - (M) hold (3.30) 



where the other periodic solution corresponding to H 2 bifurcates to the 2D motion with 
frequency $1. 

Finally, a 3D dynamical behavior can be even predicted when the following line ^5: 

L$M® (N® - M[ X) ) + 4 2 M (1) " M[ 2) )] c 2 + (3.31) 

L^M\ 2) (NP - M^) + L?hl\ l) (N^ - M^)] d 2 = 



12 lvl i y iy i ±y± \ J 12 1 y\ ^1 

is located between the lines ^3 and I4 as follows from conditions (13.25f )- f|3.27|) . 

In the post-double-Hopf bifurcation regime, the solutions exist on a torus in phase-space. 
The possible routes to chaos would most likely one of the quasiperiodic routes, i.e., either 
torus doubling via period doubling of the torus, or gradual torus breakdown, or the Ruelle- 
Takens route into chaos [15J. If the parameters are such that the system is locally very 
strongly dissipative or volume-contracting, intermittency is a possibility following bifurcations 
of the torus since the repulsion from the bifurcated (and now-unstable) torus may combine 
with volume contraction to set up a repulsion-reinjection intermittency scenario. However, 
this is less likely for most parameter sets than the previous three quasiperiodic routes. In 
exceptional cases, crises may also be possible, but again, this is much less likely for most 
system parameters than the quasiperiodic routes. 

The modified Chen system (12. 3D is next integrated for various parameters sets, showing the 
following three main behaviors: 

1. the trajectories fly off to infinity; 

2. the trajectories evolve towards a limit cycle; 

3. the trajectories evolve to a strange attractor. 

In all the numerical simulations the initial conditions are chosen to be close to the fixed 
point. The case 1., shown in Fig. [21 is realized at the parameter values a = 1,6 = 7 and 
k = 9 for which the double Hopf conditions are satisfied at c* = 8 and d* = —35. Choosing 
C2 = d,2 = 0.1 the normal form (|3.9|) - f|3. lQj) admits no stable fixed point (in particular the 
points H2 and T exist unstable, instead H\ is complex) and it is not dissipative, therefore 
the trajectories fly off to infinity. 

The case 2. occurs for the parameter choice a = 1,6 = 3, k = —6.5 (the critical parameter 
values are c* = 4 and d* = —38.5833). Choosing the second order deviation C2 = —19 and 
d 2 = —20 the system is strongly dissipative and the equilibrium point Hi exists stable. The 
other two equilibria H 2 and T are complex. Therefore our analysis predicts that the system 
states evolve towards a periodic orbit, in agreement with the simulation shown in Fig. [3[ 
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Figure 2: The time behavior of the state x for the system (|2.3[) . 



(a) 

Figure 3: The trajectories evolve towards a limit cycle. 

Finally, for the choice a = l,b = 3,k = —7 and the second order deviations C2 = — 3 
and d,2 = —3.6667, no stable equilibria exist for the normal form f)3.9[) - (|3. 10|) (in particular 
Hi exists unstable and the equilibrium points H2 and T are complex). Due to the strong 
dissipativity, the numerical simulation in Fig. [Ushows that the trajectories evolve to a strange 
attractor. 

Note that other possible route to chaos can be investigated, e.g. in [6] the post-bifurcation 
dynamics in the context of two intermittent routes to chaos (routes following either subcritical 
or supercritical Hopf or double Hopf bifurcation) are observed. 

4 Adaptive synchronization of the hyperchaotic system 

The synchronization of two chaotic/hyperchaotic systems consists into designing a controller 
or forcing in such a way that the motion of each system can be adjusted to a common dynam- 
ics. The intrinsic nature of the chaotic systems does not obey to this idea of synchronization 
due to the sensitivity to initial conditions; in fact the trajectories of two identical chaotic sys- 
tems evolve to completely different dynamical behavior when starting from different initial 
conditions (however they are close). Nevertheless, in the last decades, various synchroniza- 
tion methods have been proposed which show how it is possible to synchronize this kind of 
systems [H [281 S OS 13 M ■ 
In our case we want to realize a complete synchronization between two identical modified 




S 





(a) 



Figure 4: The strange attractor in the xyz-sp&ce. 



Chen system. The first one is the following master or driver system: 

O'iUrn %m) 



Um 



—dx r 

^rnVn 
- Xm. 



XZ m -\- Cy m Wm 

bz m k(x m + Um, 



k 2 

b 



(4.1) 



and the second one is the slave or response system: 

x s = ai(y s - x s ) + m 
ij s = -dix s - xz s + ciy. 



z s 



x s y s - b x z s 

- X s + U4 



kf 

- x s + k\z s 

01 



U2 



(4.2) 



h{x s + y s ) + u 3 



whose evolution is guided by the controllers u\,U2, u% and U4, and the parameters a\ , b\, c\ , d\ 
and k\ need to be estimated in such a way that the two systems (|4.ip and (|4.2p can be 
synchronized. 

Our goal is to design suitable control functions Ui and to find proper update rules for the 
parameters a\,b\, ci,d\ and k\ such that the response system globally synchronizes the driver 
system. This is equivalent to require that the error dynamical system, obtained by subtracting 
equations (|4.1|) from (|4.2|) : 



(•■2 



(34 



ai(y s - x s ) - a(y m - x m ) + ui 
-d\x s + dx m - xz s + x m z m + ay s 
k\ k 2 



cy ri 



W s + W T . 



, ~x s ~\~ , x m + k\z, 
b\ b 



XsUs x m y m 
x s x m U4. 

,e 2 = y s - y m , e3 



kz m + u 2 
biz s + bz m - ki(x s 



(4.3) 



where e\ = x s — 
origin, i.e. lim^oo \ \e(t)\\ = 0, e = [e x e 2 e 3 e 4 ] 



,e 4 



w. 



y s ) + k(x m + y m ) + u 3 

w m , is asymptotically stable to the 
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' e a 


= -{Vm ~ 


e-b 




< e c 


= -y m e-2 


ed 

V 


= x m e 2 
k 

= T x m e 2 




Theorem: For any initial conditions, the two systems (|4,ip and (|4.2p are globally asymp- 
totically synchronized by the following control law: 

u\ = -(hi - ai)ei - (01 — d\ — z s )e 2 - e 4 

( k k\ \ 

u 2 = -(h<2 + ci)e 2 - eie 3 - fei ( ~x m - —x s ) ^ ^ 

■"3 = -(^3 - bi)e 3 + fciej - y m ei 
u 4 = — /i 4 e 4 + e 2 

where hi are positive scalars (called control gains) and by the parameter update rules: 

Xm)e\ 

(4.5) 

- z m e 2 + (x m + V 
where e a = a\ — a, e& = b\ — b, e c = c\ — c, e<j = d\ — d, = k\ — k. 

Proof: The proof of the Theorem is based on the Lyapunov stability theory. Let us choose 
the following Lyapunov function: 

V = j (e? + el + 4 + 4 + 4 + 4 + 4 + e\ + eg) , (4.6) 

which is positive definite. Calculating the time derivative of the Lyapunov function (j4.6j) 
along the trajectories of the error system (|4.3p one obtains: 

F = eiei + e 2 e 2 + e 3 e 3 + e 4 e 4 + e a e a + + e c e c + + e^e*, 

= ei (ai(j/ s - x s ) - a(y m - x m ) + u\) + e 2 {-dix s + dx m - xz s + x m z m 

k 2 k 2 \ 
+ ciy s - cy m -w s + w m - -±x s + -^x m + k\z s - kz m + u 2 \ (4.7) 

+ e 3 (x s y s - x m y m - b 1 z s + 6z m - k\{x s + y s ) + fc(x m + y m ) + u 3 ) 
+ e 4 (x s - x m + u 4 ) + e a (ai - a) + e 6 (&i - 6) + e c (ci - c) + e d (di - d) 
+ e k (h-k). 

Substituting equations (|4.4|) and (|4.5p into the expression (j4.7j) for V and noting that the 
following equalities hold: 

Xs^s ~\~ X m Z m — Z s e± ImQi ( / L8) 

XsUs x m y m — x s e 2 + j/^6i, 

one gets: 

4 

V = — ^2 hie 2 = — e T Pe, where P = diag{^i, h 2 , h 3 , h&}. (4.9) 

i=i 
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Figure 5: Chaos synchronization: the trajectories of the error dynamical system asymptoti- 
cally converge to the origin. 

Since V is negative semidefmite ei, e2, 63, e^, e a , e^, e c , ej, e& G and from the error system 
(|4.3p it follows that e±, £3, £4 G £qo. Given A m i n (P) the minimum eigenvalue of the matrix 
P, one gets: 

t rt <-t 

A min (P)||e|| 2 dt < / e T Pedt< / -Vdt = V(0) - V(t) < V(0). (4.10) 
J Jo 

Therefore e±, e%, e^, G ^2 and the hypotheses of the Barbalat's lemma (see [20] for details) 
are satisfied. Thus lim^oo ||e(t)|| = and the proof is completed. 

□ 

To test the effectiveness of the proposed adaptive synchronization scheme, we show two 
numerical examples in which the modified Chen system has been chosen both in its chaotic 
and hyperchaotic regime. By choosing the parameters a = 1,6 = 3, c = l,d = — 49, k = —7 
the master system evolves to the chaotic attractor shown in FigJH (we have checked that only 
one Lyapunov exponent is positive). Let us assume that the starting points for the master 
and the slave systems are respectively (x m (0), y m (0), z m (0), w m (0)) = (0.1,0.3,0.01,-2) and 
M0), jfc(0), z a (p),w,(0)) = (-0.1, 0.3, -0.01, -0.2). 

Given the initial conditions for the uncertain parameters of the response system ai(0) = 
3,6i(0) = 5, ci(0) = 2, di(Q) = — 30, &i(0) = —5, the adaptive control laws (|4.4|) with control 
gains (h\, /12, ^3, ^4) = (5, 7, 6, 5) realize the chaotic synchronization of the systems (|4.ip and 
(|4.2j) as shown in Figures [5EJ 

In the second numerical test the master systems has been chosen in its hyperchaotic regime, 
with the parameters a = 36,6 = 3,c = 28, d = —16, A; = 1 (the two positive Lyapunov 
exponents are shown in Fig(T|). Let us choose the initial conditions of the master and 
the slave systems respectively as (x m (0),y m (0), z m (0) , w m (0)) = (0.3,0.03,0.1,0.2) and 
(x s (0),y s (0), z s (0), if s (0)) = (—0.1, 0.3, 0.01, —0.3) . Given the initial estimated parame- 
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Figure 6: The master and the slave system have been synchronized to the chaotic attractor 
shown in the last subfigure. The parameters of the slave system asymptotically converge to 
the master system parameter. 

ters oi(0) = 30,6i(0) = 5, ci(0) = 25,di(0) = -20,&i(0) = -1 and the control gains 
{h\,h2,h$,hi) = (5,4,6,5) the effectiveness of the hyperchaotic synchronization is shown in 
Figures EH 

5 Conclusions 

In order to understand the onset of hyperchaotic behavior recently observed in many dynam- 
ical systems, we have constructed and analyzed the generalized double Hopf normal form in 
the modified Chen system which reveals possible regimes of periodic solutions, two-period 
tori, and three-period tori in parameter space. Numerical simulations are provided showing 
agreement with the predictions from the normal form. In a future work other possible routes 
to chaos will be investigated, as further bifurcations of the post-supercritical-Hopf two- and 
three-tori via either torus doubling or breakdown. Moreover, an adaptive synchronization of 
the chaotic/hyperchaotic system trajectories has been globally realized and the effectiveness 
of this control strategy has been numerically illustrated. 
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Figure 7: Hyperchaotic synchronization: the error dynamical system asymptotically con- 
verges to the origin. 




Figure 8: The estimates of the slave system parameters converge to the master system pa- 
rameters. Both the systems (|4.ip and (|4,2p evolve to the hyperchaotic attr actor shown in the 
last sub-figure. 
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